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SECTION  I 


INTRODUCTION 

The  development  of  miniaturized  solid  state  sensors  and  associated  electronics 
has  permitted  design  of  radar  proximity  fuzes  for  use  in  small  caliber,  high 
explosive  projectiles.  However,  the  development  of  miniaturized  power  supplies 
has  not  been  commensurate  with  advances  in  micro-electronic  technology. 
Miniaturized  thermoelectric  generators  for  converting  projectile  aerodynamic 
heating  into  electricity  to  power  electronic  microcircuits  may  provide  advan  t r , 

especially  shelf-life  and  safety,  in  armament  fuzing  applications.  While  the 
results  of  this  effort  are  intended  to  be  applied  directly  to  a 20  mm  proje-'i 
fired  from  an  aircraft,  indications  are  that  there  are  many  other  applications 
in  munitions. 

1.  REVIEW  OF  PREVIOUS  ZE  DEVELOPMENT  EFFORTS 


The  Air  Force  demonstrated  the  feasibility  of  producing  a radar  proximity  fuzi 
for  a small  caliber  (20  mm)  projectile  under  Contract  F08635-71-C-0154.  As 
result  of  this  effort,  the  micro-electronics  for  a small  cw  doppler  radar 
proximity  fuze  were  developed,  and  38  prototype  fuzes  were  produced  and  test 
fired.  In  order  to  demonstrate  the  operation  of  the  electronics  circuit, 
stored  energy  in  the  form  of  a small  spin-activated  battery  was  used  as  a 
power  source. 

The  Air  Force  also  sponsored  work  under  Contract  F08635-72-C-0012  to  study, 
refine,  and  characterize  a power  converter  nose  cone  using  thermoelectric 
energy  conversion  techniques  to  transform  aerodynamic  heating  into  the  required 
electrical  power.  ^ This  contract  was  directed  primarily  towards  the  identi- 
fication of  potential  applications  for  thermoelectric  generators,  the  study  of 
thermoelectric  material  characteristics,  and  possible  manufacturing  technology 
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for  thermoelectric  generators.  It  was  apparent  from  this  preliminary  study 
that  further  research  was  needed  to  define  projectile  surface  heating  rates, 
establish  an  optimum  thermopile  design,  and  to  evaluate  mechanical  suitability 
and  projected  costs  of  such  a power  source. 

2.  OUTLINE  OF  RESEARCH  APPROACH 

The  overall  objective  of  this  Phase  I study  was  to  evaluate  the  feasibility 
of  aerodynamically  heated  thermoelectric-power  supplies  suitable  for  use  in 
small  caliber  fuzes.  To  this  end,  Battelle  conducted  a program  consisting  of 
eight  tasks  as  summarized  below 

• In  the  first  period  of  Phase  I,  the  principal  technology 
groups  associated  with  this  project  were  integrated  into 

a team  to  review  the  possible  candidate  thermopile  designs 
and  to  establish  realistic  "bounds"  for  the  subsequent 
analytical  phase.  (This  is  represented  schematically  as  Tasks 
I and  II  in  Figure  1.) 

• In  the  second  period  of  Phase  I the  thermopile  analytical 
effort  was  partioned  into  distinct  (yet  interacting)  areas 
of  consideration  such  as  aerodynamic  heating,  thermoelectric 
design  optimization,  mechanical  suitability,  cost  effectiveness, 
and  safety  (see  Tasks  III  through  VII  in  Figure  1). 

• In  the  second  phase  of  this  study  (if  continued)  BCL  plans  to 
focus  its  effort  on  prototype  hardware  development  based  on 
the  collective  results  of  the  first  phase. 

Although  the  analysis  in  the  Fhase  I study  was  partitioned  according  to  the 
principal  area  of  consideration,  each  analysis  remained  interactive  with 
respect  to  the  other  system  design  considerations,  m this  way,  the  develop- 
ment of  the  "optimum"  thermopile  progressed  through  a h'erarchy  of  design 
considerations  in  which  the  thermopile  design  becomes  further  refined  at  each 
successive  level.  As  indicated  in  Figure  1 (see  Task  III),  the  initial  level 
of  the  hierarchy  involved  that  phase  of  the  design  optimization  which  is 


2 


dictated  by  the  given  mission  and,  hence,  least  flexible.  The  next  level  of 
the  hierarchy  (see  Figure  1,  Task  IV)  involved  the  thermoelectric  design 
optimization  analysis  which  is  also  relatively  Inflexible  due  to  the  limited 
number  of  available  materials  which  meet  the  thermoelectric  energy  conversion 
requirements. 

The  discussion  which  follows  has  been  divided  into  sections  according  to  the 
tasks  outlined  in  Figure  1. 
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SECTION  II 


POWER  REQUIREMENTS  OF  PROXIMITY  FUZE 
CIRCUIT  INVOLVING  THERMOELECTRIC  POWER  SOURCE 

The  objective  of  this  task  was  to  determine  the  power  requirements  of  the 
proximity  fuze  circuit  involving  an  aerodynamlcally  heated  thermoelectric 
power  source. 

A fuze  circuit  was  developed  under  a previous  contract  involving^  power  require, 
ments  of  about  300  mw.  However,  the  proximity  fuze  circuit  power  requirements 
need  to  be  reduced  to  60  mw  or  less  in  order  to  be  feasible  for  operation  with 
an  aerodynamically  heated  thermoelectric  power  source. 

1.  CHARACTERIZATION  OF  ORIGINAL 
FUZE  CIRCUIT  DESIGN 

The  characterization  of  the  original  fuze  circuit  was  initiated  after  recu-jp: 

(2) 

of  a unclassified  version  of  the  contractor  report  entitled  "Air  to  Sur  « 
Proximity  Fuze  for  20mm  Projectile",  which  described  the  fuze  circuit  develop 
ment  and  testing. 

The  fuze  circuit  developed  in  the  previous  work  is  shown  in  Figure  2.  It 
consists  of  an  oscillator  and  an  antenna,  signal  detector,  amplifier,  impact 
switch,  a time  delay  circuit,  and  a voltage  regulator.  The  oscillator  consists 
of  transistor,  Q^,  coll,  L^,  and  the  collector/emitter  capacitance  of  transist  r 
Q , The  oscillator  output  is  coupled  to  the  antenna  and  transmitted.  The 
return  signal  is  received  by  the  quarter  wave  antenna  across  the  emitter 
resistor,  , . This  voltage  is  coupled  through  a 47  nF  capacitor  to  the 
inverting  (-)  input  of  the  operational  amplifier  while  the  non-inverting  (+) 
input  is  connected  to  the  base  of  the  transistor,  Q^,  Therefore,  this  input 
always  sees  a positive  voltage  which  never  exceeds  2.1  volts,  l.e.,  the  sum 
of  the  voltage  drops  for  the  three  diodes,  DS1.  This  permits  the  operational 
amplifier  to  act  as  a voltage  follower  and  to  operate  from  a single  battery 
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FIGURE  2.  SCHEMATIC  DIAGRAM  OF  RF  PROXIMITY  FUZE 


supply.  Thus,  the  voltage  of  sufficient  magnitude  is  received  at  the  inverting 
input  and  is  amplified;  the  output  exceeds  2.1  V (which  corresponds  to  the  sum 
of  the  voltage  drops  across  the  three  diodes,  DS2)  and  turns  on  the  silicon 
controlled  rectifier  (SCR),  which  permits  current  through  the  detonator. 

The  SCR  can  also  be  turned  on  with  the  closing  of  the  impact  switch,  SI. 

The  voltage  regulator  consists  of  one  transistor,  Qg,  its  biasing  resistors 
R10  and  R^,  and  the  Zener  diode,  VRl.  Resistor  and  capacitor  C^  were  used 
for  detection  of  detonator  firing  during  tests  and  are  not  a part  of  the  fuze 
circuit.  The  voltage  regulator  maintains  a constant  operating  voltage  for  the 
circuit.  The  time  delay  portion  of  the  circuit  consists  of  Rg,  the  fusible 
resistor,  transistor  C^,  and  resistor  R^.  A time  delay  was  needed  previously 
since  the  battery  became  active  immediately  after  shell  firing  and  the  fuze 
circuit  could  fire  the  detonator  in  the  vicinity  of  its  point  of • origin.  The 
time  delay  served  to  keep  the  fuze  circuit  in  an  unarmed  state.  This  was 
performed  through  the  use  of  Rg  and  Q^.  Current  in  Rg  turned  on 
(?2  and  kept  it  in  a saturated  state,  thus  holding  the  voltage  to  the  gate  of 
the  SCR  near  zero.  The  fusible  resistor  burned  out  after  approximately 
50  msec  and  required  300mw.  After  Rg  open  circuited,  turned  off,  and  the 
SCR  could  be  fired  by  the  operational  amplifier. 

(3) 

It  was  determined  that  the  voltage  output  of  the  voltage  regulator  was 
approximately  10  volts,  the  oscillator  frequency  was  in  the  region  of  500  Hhz> 
and  the  minimum  oscillator  power  output  was  l-3mw.  The  detonator  firing 
requirements  were  obtained  from  the  manufacturer  and  are  plotted  in  Figure  5. 

Since  this  circuit  had  been  designed  for  operation  with  a battery,  several 
changes  were  needed  to  make  it  operate  with  a thermoelectric  power  source. 

2.  DESCRIPTION  OF  MODIFIED  PROXIMITY 
FUZE  CIRCUIT 

The  proximity  fuze  circuit  needed  to  be  modified  to  have  lower  power  require- 
ments and  yet  meet  the  specified  performance  requirements.  The  first  item  to 
be  removed  was  the  fusiable  resistor.  In  the  previous  design,  the  fuzible 
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FIGURE  3.  DETONATOR  FIRING  REQUIREMENTS 


resistor  provided  a delay  in  the  arming  of  the  proximity  fuze  to  provide  a 
safe  arming  distance.  In  the  case  of  the  thermoelectric  powered  fuze,  the 
time  delay  would  occur  naturally  due  to  the  aerodynamic  heating  effects  and 
the  thermal  inertia  of  the  thermoelectric  elements.  Also,  a storage  capacitor 
was  added  to  supply  the  current  necessary  to  activate  the  detonator. 

The  minimum  voltage  for  the  circuit  is  6 volts.  This  lower  limit  is  based  on 
the  threshold  voltage  of  the  operational  amplifier.  Thus  an  operating  voltage 
of  approximately  7.5  volts  was  chosen  for  this  study.  The  power  requirements 
at  this  voltage  would  be  in  the  range  of  40  to  50mw.  Based  on  these  considerations, 
the  proximity  fuze  circuit  was  revised  and  is  shown  in  Figure  4. 

The  voltage  regulator  limits  the  applied  voltage  to  7.5  volts.  The  oscillator 
circuit  was  not  modified  since  it  had  been  optimized  in  the  earlier  study. 

Diodes,  D^,  were  added  to  provide  the  transient  to  insure  oscillator  starting 
since  the  thermoelectric  generator  is  a slowly  rising  voltage  source. 

It  must  be  remembered  that  the  revised  circuit  is  a preliminary  design  and 
has  been  neither  built  nor  tested.  Therefore,  the  numbers  given  for  power 
requirements  are  best  estimates  based  on  present  available  information. 

One  could  expect  some  changes  in  design  in  going  from  such  a theoretical  paper 
design  to  actual  hardware  and  as  the  thermoelectric  power  source  becomes  better 
characterized. 

3.  FUZE  CIRCUIT  POWER  REQUIREMENTS 

The  original  circuit,  Figure  2,  was  designed  to  operate  at  10  volts.  The 
operational  amplifier  used  is  a low  power  integrated  circuit  and  consumes  4mw; 
the  oscillator,  130mw,  and  the  voltage  regulator,  96mw,  giving  a total  of  230mw. 

No  energy  storing  capacitor  was  needed  to  fire  the  detonator  because  the  power 
supply  was  capable  of  delivering  the  necessary  current  (greater  than  lOma) . Thus., 
this  circuit  at  its  operating  voltage  point  is  not  compatible  with  a thermo- 
electric power  source. 

The  revised  circuit  shown  in  Figure  4 was  analyzed  to  determine  the  power 
requirements  at  7.5  volts,  the  lowest  voltage  at  which  the  circuit  will 

operate.  The  analysis  was  also  performed  at  a nominal  voltage  of  11  volts  (see  Figure  5). 
Analysis  of  this  circuit  was  based  on  the  requirement  for  an  oscillator  power 
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FIGURE  4.  REVISED  FUZE  CIRCUIT  (SUPPLY  VOLTAGE,  7V) 
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FIGURE  5.  FUSE  CIRCUIT  (SUPPLY  VOLTAGE,  II  V) 


output  of  l-3mw  and  on  the  aaaumptlona  of  the  efficiencies  of  the  oscillator 
and  the  coupling  to  the  transmitter  antenna.  The  oscillator  circuitry  was  not 
altered  because  of  the  previous  effort^)  expended  in  optimizing  the  oscillator 
design.  The  voltage  regulator  was  retained  because  the  thermoelectric  power 
supply  will  characteristically  change  with  elapsed  time  since  firing. 

The  study  revealed  that  the  lowest  voltage  at  which  the  fuze  circuit  will 
operate  is  six  volts.  This  lower  limit  is  due  to  the  voltage  supply  required 
by  the  operational  amplifier.  Due  to  tolerance  variations  in  Zener  diodes  and 
the  voltage  drop  across  <33,  a minimum  supply  voltage  of  7.5  volts  will  be 
needed.  The  power  requirements  at  this  lowest  voltage  are;  oscillator,  35mw; 
operational  amplifier,  2mw;  voltage  regulator,  7n w;  giving  a total  of  45mw. 

The  oscillator  is  estimated  to  generate  a signal  with  a power  of  4~6mw.  Since 
neither  the  frequency  nor  oscillator  efficiency  was  available,  this  approximation 
becomes  necessary. 


12 


SECTION  III 


IDENTIFICATION  OF  CANDIDATE  THERMOPILE  CONFIGURATION 
AND  ASSOCIATED  MATERIALS 

Based  on  the  fuze  circuit  power  requirements  identified  in  the  previous  section 
of  this  report,  e preliminary  thermoelectric  analysis  was  performed.  The  thermo- 
electric analysis  was  performed  parametrically  with  respect  to  choice  of  thermo- 
electric material,  length  of  thermoelectric  elements,  and  hot- junction  temperature. 

The  results  of  the  preliminary  thermoelectric  analyses  are  sunmarized  in  Table  I 
and  Indicate 

• Number  of  thermoelectric  couples  required  for  selected 
temperature  difference  across  thermopile 

• Cross-sectional  area  of  thermoelements  corresponding  to 
selected  temperature  differences  and  thermopile 
configuration 

• Number  of  thermopile  units  or  modules  required  to  meet 
fuze  voltage  and  current  requirements. 

1.  SELECTION  OF  RANGE  FOR  PARAMETERS 

The  selection  of  each  of  the  parameters  summarized  in  Table  I was  based  on  the  estimated 
practical  limits  of  materials  properties,  fabricabillty , and  dimensional  control. 

The  rationale  underlying  the  selection  of  the  range  for  each  parameter  is  briefly 
summarized  below. 

a.  Thermoelectric  Materials 

There  are  five  general  classes  of  materials  that  exhibit  thermoelectric 
characteristics:  metals,  Scmimetals,  semiconductors,  semi insuls tors , and  insulators. 

The  metals  are  too  low  in  Seebeck  coefficient  while  too  high  in  thermal  conductivity 
to  be  of  practical  significance.  In  contrast  the  semiinsulators  and  Insulators, 
while  high  in  Seebeck  coefficient,  are  too  high  in  electrical  resistivity.  Of 
these  five  classes  of  materials,  only  the  semiconductors  and  semimetals  are  suitable 
for  the  present  energy  conversion  application.  Within  the  semiconductor  class, 
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TABLE  1. 

SUMMARY  OF  THERMOELECTRIC  CONVERTER  CHARACTERISTICS  FOR 
SELECTED  SEMICONDUCTOR  AND  SEMIMETAL  MATERIALS 
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those  materials  offering  the  highest  energy  conversion  efficiency  include  bismuth 
telluride,  lead  telluride,  and  silicon  germanium.  Chromel-P/Alumel  and  Chromel-P/ 
Constantan  combinations  are  representative  of  the  semlmetal  class  of  materials. 

b.  Temperature  Difference 

The  selection  of  the  range  of  temperature  differences  between  the  hot- 
junction  and  cold- junction  of  the  thermopile  was  based  on  the  minimum  required 
temperature  difference  for  practical-sized  thermopiles  and  the  maximum  anticipated 
temperature  difference  based  on  the  available  aerodynamic  heating  rates. 

c.  Thermoelement  Length 

The  upper  limit  of  the  range  selected  for  this  parameter  was  determined 
by  the  overall  dimensions  of  the  projectile.  The  lower  limit  on  the  thermoelement 
length  is  influenced  by  the  practical  limits  of  the  proportionately  smaller  cross- 
sectional  area  of  the  thermoelements  as  well  as  the  increasing  proportion  of  power 
dissipated  at  the  thermoelement/electrode  junctions.  In  addition,  at  decreasing 
thermoelement  lengths,  hence  cross-sectional  areas,  the  thermal  heat  flux  require- 
ments will  exceed  the  available  aerodynamic  heat  flux. 

d.  Thickness  of  Thermoelements 

Based  on  the  preliminary  thermoelectric  analyses  using  the  averaged- 
(4) 

parameter  technique,  it  was  apparent  that  the  cross-sectional  dimensions  of 
the  individual  thermoelements  (for  all  materials  considered)  were  beyond  the 
limits  of  cost-effective  fabrication  by  the  "discrete  element"  approach.  One 
alternative  fabrication  method  involves  the  vapor  deposition  of  a thin  film  of 
the  thermoelement  on  an  appropriate  substrate.  For  the  purposes  of  the  present 
comparative  analysis,  film  thicknesses  of  0.00254  cm  (0.001  in.)  and  0.00508  cm 
(0.002  in.)  have  been  selected. 

2.  SELECTION  OF  THERMOPILE  CONFIGURATION 

With  the  range  for  thermopile  design  parameters  established,  the  next  step  was 
to  Identify  a feasible  thermopile  configuration.  In  order  to  provide  the  required 
substrate  surface  area,  a concept  featuring  the  use  of  a multiplicity  of  disc- 
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shaped  modules  (see  Figure  6)  has  been  evolved.  In  this  scheme,  a series 
arrangement  of  thermoelements  is  deposited  on  both  sides  of  the  disc-shaped 
substrate.  This  overall  module  design  scheme  was  used  in  the  comparative  thermo- 
electric analyses  discussed  next.  The  physical  location  of  the  disc  modules  in 
the  projectile  tip  is  shown  in  Figure  7. 

3.  SELECTION  OF  CANDIDATE  THERMOPILE  DESIGN 

The  results  of  the  comparative  thermoelectric  analysis  summarized  in  Table  1 
indicate  that  the  bismuth  telluride  thermopile  features  a number  of  advantages 
including 

• Fewest  number  of  couples  required 

• Fewest  number  of  disc-shaped  modules  required 

• Lowest  heat  flux  requirement 

• Practical  element  width  from  standpoint  of  physical  masking. 

The  semimetals  are  attractive  from  the  standpoint  of  a simpler  materials  system 
(no  dopant  required)  and  a lower  base  cost  for  the  materials.  However,  the 
thermopiles  involving  semimetals  involve  very  small  thermoelement  widths,  thus 
complicating  the  vapor  deposition  process  and  associated  masking  requirements. 

This  analysis  indicates  that  bismuth  telluride  is  the  most  promising  thermoelectric 
material  for  the  present  20-mm  projectile  power  supply  requirements.  Based  on 
these  results,  more  detailed  analysis  of  the  overall  thermopile  performance  was 
performed  and  is  reported  in  Section  V. 
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FIGURE  7.  CROSS-SECTIONAL  VIEW  OF  20  mm 
PROJECTILE  INVOLVING  DISC 
THERMOELECTRIC  MODULES 
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SECTION  IV 


AERODYNAMIC  HEATING  ANALYSIS 

The  objective  of  the  aerodynamic  heating  analyala  la  to  establlah  the  maximum 
and  minimum  aurface  heating  ratea  experienced  by  a typical  20  mm  projectile. 
Aircraft  altitude  and  velocity  are  to  be  Included  aa  variables,  and  the  heating 
ratea  are  to  be  determined  aa  a function  of  poaitlon  along  the  projectile 
surface.  The  results  of  the  analysis  will  be  used  as  Input  to  the  thermopile 
design  optimization. 

1.  BLUNT  CONE  CONFIGURATION 

The  blunt  cone  configuration,  typical  of  a 20  mm  projectile  shown  in  Figure  8, 
produces  a rather  complex  flow  field.  The  curvature  of  the  bow  shock  wave 
generated  by  the  blunt  nose  of  the  projectile  causes  the  entropy  to  vary  from 
streamline  to  streamline.  Due  to  the  growth  of  f ie  boundary  layer  across 
streamlines,  the  fluid  properties  in  the  boundary  layer  vary  with  surface 
distance.  This  effect  is  propagated  far  back  on  the  conical  portion  of  the 
body.  As  a result,  boundary  layer  heating  cannot  be  accurately  predicted  by 
simple  analytical  techniques  which  rely,  for  example,  on  the  principle  of 
boundary  layer  similarity. 

After  examining  the  various  available  computer  programs  with  the  capability 
to  calculate  convective  heating  to  blunt  cones,  the  NOL  RETAS  code  was  chosen. 
This  program  has  been  developed  over  several  years  by  personnel  at  the  Naval 
Ordnance  Laboratory to  compute  drag  and  convective  heating  on  spherically 
blunted  cones.  The  pro yic tile  trajectory  is  computed,  and  drag  and  distributed 
convective  heating  are  c ..lculated  at  selected  points  along  the  trajectory. 

The  code  is  ideally  suiced  to  this  subtask. 

2.  LAMINAR  BOUNDARY  LAYER  HEAT  TRANSFER 

The  NOL  RETAS  Code  utilises  the  energy- integral  approach  to  determine  the 
Stanton  number  distribution  along  the  surface  of  a blunted  cone  at  zero  angle- 
of-attack.  The  energy-integral  approximation  is  similar  to  the  Von  Karman 
momentum- integral  approach.  An  energy-integral  equation  is  developed  by 
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FIGURE  8.  SCHEMATIC  OF  BLUNT  CONE  FLOW  FIELD  TYPICAL  OF  20nm  PROJECTILE 


Integrating  the  energy  equation  through  the  boundary  layer.  The  reaultlng 
equation  la  then  solved  by  using  a reference  temperature  method  and  an  lncompress lb le, 
flat-plate  Stanton  number  relation.  The  solution  assumes  perfect  gas  flow  over 
the  cone.  Complete  details  of  this  method  are  given  In  References  5 and  6. 

Once  the  boundary  layer  properties  are  determined,  the  local  skin  friction 
coefficient  and  the  heat-transfer  rate  are  found.  The  skin  friction  coefficient 
Is  Integrated  over  the  body  and  contributes  to  the  total  drag  calculation. 

The  local  heat  transfer  rate  is  given  by 

q„  * P,U,C  (T,*-T  ) St 
K1  1 p'  1 w' 

where  St  Is  the  Stanton  number,  p.  and  U.  are  the  boundary  layer  edge  density 

^ ^ * 

and  velocity,  respectively,  Ty  is  the  wall  temperature,  and  T^  Is  the  reference 
temperature. 

The  surface  heat  transfer  rates  calculated  in  this  manner  have  been  compared 
with  the  blunt  cone  measurements  of  Reference  7.  The  results  are  presented  in 
Figure  9 as 

\ VK(Ti*-V 

R. ' R. 

versus  */rn.  The  calculations  compare  very  well  with  the  data  on  the  conical 
part  of  the  body. 

3.  PROJECTILE  SPIN  EFFECTS 

The  high  rate  of  spin  of  a typical  20  mm  projectile  produces  magnus  forces  and 
moments  which  are  Important  to  the  flight  dynamics  of  the  projectile.  It  Is 
also  possible  that  the  spinning  of  a projectile  will  influence  the  character 
of  the  boundary  layer  so  that  boundary  layer  heating  will  be  affected. 

The  first  order  effect  of  spin  can  be  estimated  by  comparing  the  true  surface 
velocity.  Including  spin,  with  the  projectil*.  axial  velocity.  The  true  surface 
velocity  Is  obtained  by  vectorially  adding  the  axial  and  tangential  velocities 
of  the  surface.  For  an  axial  velocity  (muzzle  velocity)  of  4000  ft/sec  and  a 
projectile  spin  of  120,000  rpm,  the  true  surface  velocity  is 


* s/rn  refers  to  distance  along  projectile  surface (s)  from  the  stagnation  point 
divided  by  the  nose  radius,  rn. 
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FIGURE  9.  COMPARISON  OF  EXPERIMENTAL  DATA  WITH  THEORY  FOR 
LAMINAR  FLOW  OVER  A BLUNTED  CONE,  REF  1 


2 

V * v2  . , + — * 4020  ft/aec, 

a axial  2 ’ 

where  u)  la  the  projectile  angular  velocity  and  d la  the  projectile  diameter, 
which  ia  within  1/27.  of  th*>  axial  velocity.  Thus,  as  a first  order  Influence, 
projectile  spin  can  be  neglected. 

Evidence  supporting  this  conclusion  ia  reported  in  Reference  8.  In  Reference 
8,  pitot  pressure  surveys  were  made  of  the  laminar  boundary  layer  on  an  ogive- 
cylinder  model,  at  Mach  3.02  and  4.10.  Boundary  layer  properties  were  compared 
between  zero  spin  and  a spin  rate  relative  to  axial  velocity  of  twice  that  for 
the  20  mm  projectile  discussed  above.  The  data  confirm  the  earlier  conclusion 
that  the  zero  angle-of-attack  boundary  layer  velocity  profiles  do  not  appear 
to  be  affected  significantly  by  model  spin. 

4.  BOUNDARY  LAYER  TRANSITION 

The  blunt  cone  configuration  is  generally  used  for  supersonic  and  hypersonic 
projectiles  because  it  is  an  optimum  shape  with  regard  to  low  drag  projectiles. 
Blunting  a sharp  cone  effectively  reduces  drag  in  two  ways.  First,  the  bow 
shock  wave  produced  by  the  blunt  nose  significantly  reduces  the  flow  velocity 
along  the  cone  boundary  layer  edge  below  the  value  that  would  exist  on  a sharp 
cone.  This  reduces  the  skin  friction  drag  (as  well  as  boundary  layer  heating), 
which  is  the  major  contributor  to  the  total  projectile  drag.  Secondly,  blunting 
a cone  has  been  shown  to  inhibit  boundary  layer  transition  from  laminar  to 
turbulent  flow.  A turbulent  boundary  layer  can  produce  3 or  4 times  the  skin 
friction  of  a laminar  boundary  layer. 

The  reason  nose  blunting  forstalls  transition  is  still  uncertain  and  will  not 
be  discussed  here.  There  is,  however,  a great  deal  of  boundary  layer 
transition  data  for  blunt  cones,  for  both  free  flight  and  wind  tunnel  experiments. 
Tetrevin^  has  developed  a correlation  parameter,  based  on  theoretical  consid- 
erations, which  successfully  correlates  the  results  from  48  ballistic  range  tests 
for  cones.  Tetrevinfe  correlation  predicts  transition  at  a. length  Reynolds 
number  of 
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where 


4/3 


T 2/3 

C = Y [l/2(T  /T  +1)  4.  0.0037  M2]  Re  . 

1e  9 

I'l: -j  data  and  correlation  are  presented  in  Figure  10*  The  results  are  reasonably 
well  described  by  the  correlation.  This  correlation  has  been  used  to  predict 
t e occurrence  of  transition  or  lack  of  transition  and  the  subsequent  effect 
oi.  aerodynamic  heating  for  the  20  mm  projectiles  analyzed  during  this  program. 

£ RESULTS  OF  BOUNDARY  LAYER  HEATING  CALCULATIONS 

The  range  of  parameters  over  which  the  convective  heating  to  20  mm  projectiles 
must  be  considered  is  presented  in  Figure  11.  Shown  on  the  altitude-Mach 
number  map  is  the  range  in  these  parameters  a projectile  might  experience  upon 
being  fired  from  a F-lll,  a helicoptor,  and  from  the  rear  of  a B-52.  Super- 
irr.pcsed  on  the  map  are  curves  of  constant  convective  heating  rate.  Also,  a 
curve  of  constant  Reynolds  number,  equal  to  20  x 106,  is  given.  To  the  right 
cf  this  line  it  is  possible  that  transition  to  turbulent  boundary  layer  flow  will 
occur  on  parts  of  the  projectile. 

A3  an  example  of  the  results  obtained  from  analyses  of  the  convective  heating 
re  tea  on  the  20  mm  projectile,  detailed  calculations  for  the  nominal  case  (indi- 
cated by  a * in  Figure  11)  are  presented  in  Figures  12  and  13.  A 0,20-lb 
projectile,  fired  at  3500  ft/sec  from  3000  ft  at  an  angle  of  -20  degrees  from 
e e horizontal,  was  considered  for  the  nominal  case.  The  projectile  range  and 
ulocity  history  are  shown  in  Figure  12.  In  Figure  13,  the  distributed  convective 
eating  to  the  projectile  surface  is  presented  at  three  times  during  the  first 
second  Vt  flight.  *The‘  Mat  transfer  raLe- is  seen  to  drop  of f "dramatically  in 
this  first  second  due  to  the  decrease  in  projectile  velocity.  At  some  time 
c"  o projectile  velocity  will  be  so  low  that  the  projectile  will  experience 
ceding  by  the  free  stream  instead  of  heating. 
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FIGURE  11.  ALTITUDE- MACH  NO.  MAP  FOR  20nm  PROJECTILE  FIRED  FROM  VARIOUS  VEHICLES 


Altitude,  1000  ft  Altitude,  1000  ft 


20  mm  Projectile  Heating 
Range  and  Velocity  History 
Projectile  Weight  = 0.20  lb 


FIGURE  12.  RESULTS  OF  TRAJECTORY  CALCU- 
LATION FROM  RETAS  CODE 
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Heat  Flux,  Btu/ft^  sec 


20  inn  Projectile  Heating 

Rn  - 0.125  In.;  ■ 15  degrees 


Altitude 
Muzzle  Velocity 
Firing  Angle 
Elapsed  Time 


3000  ft 
3500  ft/sec 
-20  degrees 
t 


t = 0 sec 
V = 3500  ft/sec 


t = 0.35  sec 
V = 2864  ft/sec 
position  of 
convertor 


1.25  sec 
1823  ft/sec 


Distance  along  Surface  S/Rn 


FIGURE  13.  DISTRIBUTED  HEATING  ON  20  MM  PROJECTILE 
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Calculations  (see  Figure  14)  of  heat  transfer  rates  were  performed  for  the  20-ran 
projectile  at  representative  points  throughout  the  map  of  Figure  4.  It  was 
found  that  laminur  boundary  layer  flow  prevailed  except  at  one  small  region  of 
the  F-lll  range.  Transition  occurred  near  the  base  of  the  cone  (s/rn  * 13)  at 
speeds  greater  than  Mach  4 and  altitudes  of  3000  ft  and  less.  However,  in  all 
cases  relamlnarization  occurred  as  the  projectile  slowed  down,  such  that  by 
t ■ 0.2  sec  the  flow  was  completely  laminar. 

The  aerodynamic  heating  rate  for  a 30-nra  projectile  is  included  in  Figure  14 
and  illustrates  the  effect  of  a larger  projectile  mass  on  the  velocity  decay 
rate,  hence,  'iUo  »ynamic  heating  rate.  The  performance  of  thermoelectric 
convertor  for  a 30-mm  projectile  will  be  included  in  a later  section  of  this 
report  along  with  selected  cases  involving  20-ram  projectiles. 

Substantiatin'  f the  absence  of  transition  on  the  cone  of  a 20-nrni  projectile 
is  found  in  weference  10.  In  Reference  10  the  temperature  history  of  various 
components  of  a 20-mm  projectile  was  obtained  by  transient  heating  tests  in  a 
Mach  3 blowdown  tunnel.  It  was  found  that  turbulent  flow  existed  on  the 
shell  casing,  downstream  from  the  cone,  at  a Reynolds  number  per  foot  of 
20  x 106.  Below  this  Reynolds  number,  laminar  flow  existed  on  the  entire 
shell.  Since  it  would  require  a slightly  higher  value  of  Reynolds  number  to 
move  transition  onto  the  cone,  it  is  felt  that  the  results  of  Reference  10  are 
in  agreement  with  the  transition  criteria  used  in  this  study. 
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FIGURE  14.  AERODYNAMIC  HEATING  OF  THE  PROJECTILE  AT  A POINT  1 
FROM  TIP  FOR  A FIRING  ANGLE  OF  -20  DEGREES 


SECTION  V 


OPTIMIZATION  OF  THERMOPILE  DESIGN 

The  results  of  the  thermopile  materials /design  study  discussed  in  Section  III 
were  were  translated  into  a thermal  model  for  the  purpose  of  transient  thermal 
analysis  of  the  thermopile  under  simulated  conditions  of  firing.  The  thermal 
model  for  the  disc  thermoelectric  module  configuration  is  shown  in  Figure  15. 
The  dimensions,  thermal  connections,  and  properties  associated  with  the  net- 
work of  nodes  (subregions)  shown  in  Figure  15  provide  the  input  data  for  the 
TRUMP*  heat  transfer  computer  program. 

The  input  data  for  this  thermal  analysis  include: 

• Dimensions  of  each  subregion 

• Thermal  conductivity,  specific  heat,  density  associated 
with  each  material  in  the  system 

• Heat  transfer  coefficients  associated  with  internal 
and  external  interfaces 

• View  factors  and  surface  emittances  at  internal  and 
external  interfaces 

• Surface  heat  flux  as  a function  of  time 

• Initial  temperature  of  subregions. 

The  output  data  provided  by  the  thermal  analysis  include 

• Temperature  of  each  subregion  as  a function  of  time 

• Heat  transferred  in  and  out  of  each  subregion. 

The  calculated  temperature  profile  of  the  thermoelectric  converter  as  a 
function  of  time  provides  the  basis  for  estimating  the  output  power  of  the 
thermoelectric  converter  as  a function  of  time.  These  analytical  results 
thus  provide  the  rationale  for  such  design  tradeoffs  as  increasing  (de- 
creasing) the  "thermal  storage"  capability  of  the  disc  converter  and  intro- 
ducing electrical  power  storage  for  the  latter  phase  of  the  trajectory. 


* See  Appendix  for  description  of  TRUMP  computer  program. 
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FIGURE  15.  THERMAL  MODEL  OF  DISC 

THERMOELECTRIC  CONVERTER 
CONCEPT 
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1.  DESCRIPTION  OF  THERMAL  MODEL 

The  thermal  model  for  the  thermoelectric  disc  converter  Is  Illustrated  In 
Figure  15.  The  first  step  In  the  development  of  the  thermal  model  Involved 
partitioning  the  disc  module  Into  the  smallest  sector  of  symmetry  with  respect 
to  the  temperature  profile  of  the  disc  modu’e.  In  the  present  case,  the 
smallest  sector  of  symmetry  Involved  one-half  of  the  disc  thickness  with  an 
angular  sector  which  Includes  one  thermoelectric  element  plus  the  region 
between  two  adjacent  thermoelectric  elements  (see  Figure  15).  Once  the 
smallest  sector  of  symmetry  was  Identified,  the  next  step  involved  sub- 
dividing this  sector  of  the  disc  Into  subregions.  These  subregions  were 
interconnected  by  equivalent  "thermal  resistances"  based  on  the  temperature- 
dependent  thermal  conductivity  of  the  materials  involved,  the  dimensions  of 
the  subregions,  and  the  Interface  thermal  impedances  (if  any)  between  two 
given  subregions.  The  heat  capacitance  of  each  subregion  was  obtained  from 
the  product  of  the  subregion  volume  and  the  temperature-dependent  specific 
heat  of  the  material  occupying  that  subregion. 

The  assigned  thermal  capacitances  and  thermal  coupling  between  subregions 
thus  comprise  the  thermal  model  used  in  the  present  study.  The  parameters 
for  the  thermal  analysis  of  the  connector  included 

• Substrate  material 

• Time-dependent  aerodynamic  heating  rate  (simulating 
various  altitude/muzzle  velocity  combinations) 

• Hot  strap  thickness. 

The  constants  and  property  data  used  in  this  thermal  analysis  included 

• Bismuth  Telluride  thermoelectric  element  ( k = 0.0162 
watts/cm-C,  Cp  =0.0107cal/g) 

• Steel  projectile  core  (k  = 0.512  watts/cm. C, 

Cp  = 0.460cal/g) 

• Thermal  impedance  at  the  core/disc  interlace  = 0 
(watts/cm  C) 

2 -1 

• Thermal  impedance  at  all  other  interfaces  = 0(watts/cm  . C) 

• Dimensions  as  shown  in  Figure  15. 
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In  this  model,  0.005  cm  thick  copper  hot  straps  (see  Figure  6)  were  included 
in  order  to  improve  the  transfer  of  aerodynamic  heat  from  the  outer  surface 
of  the  disc  to  the  thermoelectric  element.  As  shown  in  Figure  6,  the  copper 
hot  strap  was  positioned  over  the  outer  face  of  the  disc  and  extended  or 
"wrapped"  over  the  side  of  the  disc  and  in  thermal/electrical  contact  with 
the  thermoelectric  element. 

2.  THERMAL  ANALYSIS  OF  DISC  CONVERTER  FOR 
SELECTED  SUBSTRATE  MATERIALS 

The  next  step  in  the  optimization  of  the  thermoelectric  converter  design 
involved  identifying  preferred  substrate  (disc)  material(s).  The  criteria 
for  selecting  substrate  materials  included 

• Low  thermal  conductivity 

• Capability  of  operation  at  temperatures 
up  to  400  F (205  C) 

• High  tensile  strength. 

In  these  analyses  only  the  reference  case  was  considered,  i.e.,  20  mm 
projectile  fired  at  3000  foot  altitude  with  a muzzle  velocity  of  3500  ft/sec. 
A parametric  analysis  involving  three  substrate  materials  (alumina,  Pyroceram 
9606,  and  Kapton)  was  performed.  Ihe  results  of  this  analysis  are 

summarized  in  Table  II  and  indicate  that  the  substrate  (disc)  thermal  con- 
ductivity must  be  an  order  of  magnitude  lower  than  the  thermoelement  thermal 
conductivity  in  order  to  develop  adequate  temperature  differences  (e.g., 

50  C or  more)  across  the  thermoelements.  This  is  particularly  important 
when  one  considers  that  the  thickness  of  the  disc  is  seven  times  larger  than 
the  thermoelectric  elements. 

With  identification  that  the  disc  converter  design  involving  Kapton  allowed 
adequate  temperature  differences  to  develop  across  the  thermoelements, 
additional  parametric  analyses  were  performed  to  evaluate  temperature 
difference  across  the  converter  and  corresponding  output  power  as  a function 
of  time  for  selected  time-dependent  aerodynamic  heating  rates.  These  results 
are  discussed  next. 
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TABLE  II. 

COMPARISON  OF  PEAK  TEMPERATURE  DIFFERENCES  ACROSS 
THERMOELEMENTS  FOR  SELECTED  SUBSTRATE  MATERIALS 


Substrate 

Material 

Thermal  Conductivity,  Watts/cm-C 
Substrate  Thermoelement 

Peak 

Temperature 
Difference, K 

Alumina 

0.28 

0.016 

10 

Pyroceram  9606 

0.037 

0.016 

34 

Kapron 

0.002 

0.016 

105 
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3.  TEMPERATURE  DIFFERENCE  ACROSS  CONVERTER 
AS  A FUNCTION  OF  TIME 

A total  of  r.ine  cases  were  next  analyzed  to  determine  the  temperature  dif- 
ference developed  across  the  thermoelements  as  a function  of  time  for  selected 
altitudes  and  muzzle  velocities.  The  cases  analyzed  are  summarized  in  Figure 
11  in  terms  of  aircraft  altitude  and  projectile  Mach  number.  The  results  of 
these  analyses  are  summarized  in  Figure  16  and  indicate  that  significant 
temperature  differences  can  be  achieved  in  all  but  the  low  muzzle  velocity 
case.  This  case  corresponds  to  the  rear-firing  guns  on  a B-52  aircraft.  The 
results  were  next  used  to  estimate  the  thermoelectric  couple  output  power 
as  a function  of  time  under  matched  load  conditions. 

4.  OUTPUT  POWER  AS  A FUNCTION  OF  TIME 

Of  the  nine  cases  considered  above,  five  were  further  analyzed  to  deter- 
mine the  couple  open-circuit  voltage  and  output  power  as  functions  of  time. 
This  analysis  was  performed  as  illustrated  below  for  each  leg  of  the  thermo- 
couple : 

7 

Open  circuit  voltage,  V = £ S , (T)  6T. 

oc  i=i  i 

L 7 

Internal  resistance,  R.  = — E p.(T) 

int  A . . i 
i=i 

(V  )2 
* oc 

Output  power  , P = 

int 

The  summation  is  taken  over  the  seven  segments  which  make  up  the  overall 
thermoelement  (see  Figure  14). 

The  results  of  these  analyses  are  summarized  in  Figures  17  and  18  and  indicate 
that  an  adequate  output  power  per  couple  should  be  produced  by  aerodynamic 
heating  for  elapsed  times  of  one  second  or  more.  Note  that  at  the  higher 


* Matched  load  conditions. 
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Elapsed  Tine,  seconds 

FIGURE  16.  TEMPERATURE  DIFFERENCE  ACROSS  CONVERTOR  AS  A FUNCTION  OF  TIME 


FIGURE  17.  CONVERTOR  OPEN  CIRCUIT  VOLTAGE  AS  A FUNCTION  OF  TIME 


sa^B/vpinm  ‘aaftoj  qndqno  aidnoo 
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FIGURE  18.  CONVERTOR  OUTPUT  POWER  AS  A FUNCTION  OF  TIME 


velocity  (4000  ft/sec)  a significantly  larger  output  power  owing  to  the  high 
aerodynamic  heating  rates  is  generated.  Clearly,  for  higher  muzzle  velocities, 
there  is  some  incentive  to  produce  excess  electrical  power  during  the  early 
period  of  the  trajectory  and  store  it  for  use  in  the  later  periods  when  velo- 
city and,  hence,  output  power  have  decayed  substantially. 

Without  energy  storage,  the  number  of  couples  required  to  produce  40  milli- 
watts output  power  is  inversely  proportional  to  the  minimum  couple  output 
power  encountered  during  the  trajectory  lifetime.  This  relationship  is 
illustrated  in  Figure  18  and  indicates  that,  for  3000  feet  altitude  and  3500 
ft/sec  muzzle  velocity,  the  number  of  couples  required  for  a trajectory  life- 
time of  2.5  second  is  approximately  220.  However,  with  suitable 
energy  storage,  this  couple  requirement  might  be  reduced  to  150  couples. 

The  options  for  energy  storage  will  be  discussed  later  in  this  section. 

5.  SENSITIVITY  OF  THERMOPILE  TEMPERATURE  DIFFERENCE 
AND  OUTPUT  POWER  TO  MUZZLE  VELOCITY 

Additional  aerodynamic  heating  analyses  were  performed  to  provide  projectile 
surface  heating  rates  as  a function  of  muzzle  velocity  and  elapsed  time  after 
firing.  In  these  analyses,  the  altitude  (3000  feet)  and  firing  angle  (-20 
degrees)  were  maintained  constant  for  muzzle  velocities  ranging  from  3500  to 
4800  ft/sec.  The  purpose  of  this  limited  parametric  analysis  was  to  estimate 
the  sensitivity  of  the  aerodynamic  heating  rate  and  temperature  difference 
to  the  projectile  muzzle  velocity  over  the  period  of  the  projectile  trajectory. 
Accordingly,  the  sensitivity  of  the  related  .emperature  difference  across 
the  thermopile  can  be  translated  to  the  open  circuit  voltage  and  output  power 
as  a function  of  muzzle  velocity.  The  results  of  this  sensitivity  analysis 
thus  provided  the  basis  for  translating  20  mm  projectile  muzzle  velocity  distri- 
butions into  thermopile  output  power  and  peak  hot-junction  temperature 
distributions . 

The  results  of  this  parametric  analysis  are  summarized  in  Table  III  and  indi- 
cate that  the  temperature  difference,  AT,  and  output  power  (per  couple), 

P,  are  related  to  the  muzzle  velocity,  V,  as  follows  (after  elapsed  time  of 
0.3  second): 
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A(TH  - TC)  = 0.1565  (AV)  (3500  < V < 4800) 

AP  = 0.0012  (AV) 

P * output  power  in  milliwatts 

TC,  TH  = temperatures  in  C 

V = velocity  in  ft/sec. 

(U) 

A review  of  a reported  muzzle  velocity  uniformity  study  revealed  that  20  mm 
ammunition  based  on  M246  IED  test  firings  in  a XM163  automatic  cannon  can  be 
characterized  by  a muzzle  velocity  of  3380  + 50  ft/sec.  Based  on  the  above 
findings,  this  2.7  percent  change  in  velocity  would  correspond  to  a couple 
output  power  variation  (after  elapsed  time  of  0.3  second)  of  approximately 
22  percent. 

The  above  results  can  now  be  translated  into  operational  ’.equirements  for 
aircraft  deploying  the  subject  20  mm  projectiles.  Specifically,  the  range 
of  aircraft  speed  suitable  for  use  with  aerodynamical ly  heated  thermoelectric 
(bismuth  telluride)  power  supplies  can  be  estimated  by  assuming  a nominal 
muzzle  velocity  of  3380  ft/sec  based  on  a stationary  aircraft  cannon.  The 
results  summarized  in  Table  III  indicate  that  a minimum  effective  muzzle 
velocity  (muzzle  velocity  plus  aircraft  speed)  of  3500  ft/sec  is  required 
in  order  to  provide  adequate  output  power  up  to  at  least  a 2.5  second  trajec- 
tory (for  no  energy  storage).  At  the  other  extreme,  the  maximum  effec- 
tive muzzle  velocity  is  4500  ft/sec  to  avoid  exceeding  the  upper  operating 
temperature  limit  of  the  disc  thermoelectric  converter.  These  two  extremes 
translate  to  maximum  and  minimum  aircraft  speeds  of  1120  ft/sec  (760  mph) 
and  120  ft/sec  (80  mph),  respectively  at  a firing  altitude  of  3000  feet. 

6.  SENSITIVITY  OF  THERMOPILE  TEMPERATURE  DIFFERENCE 
AND  OUTPUT  POWER  TO  ALTITUDE  AT  FIRING 

At  higher  altitudes,  the  upper  and  lower  limit  aircraft  speeds  will  be  reduced 
due  to  the  lower  aerodynamic  heating  rates  encountered  in  the  less  dense  air 
which  lead  to  a delayed,  but  larger,  temperature  difference.  For  example, 
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TABLE  III. 

VARIATION  IN  TEMPERATURE  DIFFERENCE  AND  OUTPUT 
POWER  AS  A FUNCTION  OF  MUZZLE  VELOCITY^ 
FOR  20-MM  PROJECTILE 


Elapsed  Time 
Since  Firing, 
seconds 

Effective 

Muzzle  Velocity. 

feet/second 

3500 

3700 

4000 

It 

plcj 

AT 

P 

AT 

P 

0.3 

146 

0.339 

167 

0.511 

208 

0.752 

0.6 

167 

0.499 

194 

0.647 

236 

0.912 

0.9 

166 

0.483 

189 

0.612 

231 

0.850 

1.2 

159 

0.441 

181 

0.552 

219 

0.758 

1.5 

149 

0.384 

170 

0.482 

205 

0.665 

1.8 

139 

0.336 

156 

0.406 

187 

0.552 

2.1 

130 

0.293 

144 

0.346 

178 

0.495 

2.5 

119 

0.245 

133 

0.296 

162 

0.411 

(a)  Based  on  an  altitude  of  3000  feet  and  firing  angle 
of  -20  degrees. 


(b)  Temperature  difference  across  thermoelectric  convertor 
in  centigrade  degrees. 

(c)  Electrical  output  power  per  couple  in  milliwatts. 


42 


compare  Che  thermopile  temperature  difference  as  a function  of  time 
since  firing,  shown  in  Figure  16,  for  two  cases--15000  ft  and  3000  ft  altitude 
at  3500  ft/sec  effective  muzzle  velocity.  The  lower  altitude  (3000  ft)  case 
illustrates  the  effect  of  the  increased  aerodynamic  drag  in  the  denser  air-- 
resulting  in  an  early  peak  in  the  temperature  difference  (see  Figure  16) 
followed  by  a steep  decay  in  the  temperature  difference  during  the  remainder 
of  the  trajectory.  In  contrast,  the  high  altitude  (15000  ft)  case  reflects 
a lower  aerodynamic  drag--resulting  in  a relatively  slow  increase  in  the 
temperature  difference  across  the  thermopile  during  the  first  1.0  second  of 
the  trajectory  (see  Figure  16). 

It  is  interesting  to  note  that,  from  the  above  comparison  of  thermopile 
performance  for  low  and  high  altitudes,  the  use  of  energy  storage  (to  extend 
the  useful  range  of  the  projectile)  is  effective  only  at  lower  firing  alti- 
tudes. At  these  lower  altitudes,  excess  power  is  produced  during  the  early 
period  of  the  trajectory  whereas  at  higher  altitudes  the  major  fraction  of 
the  power  produced  occurs  later  in  the  trajectory. 

7.  CIRCUIT  ANALYSIS  OF  FUZE  CIRCUIT  POWERED 
BY  THERMOELECTRIC  CONVERTER 

To  examine  the  performance  of  the  fuze  circuit  with  the  thermoelectric  genera- 
tor, the  fuze  circuit  was  analyzed  with  the  computer  program,  SCEPTRE*. 
SCEPTRE  is  a circuit  analysis  program  capable  of  determining  the  transient 
and  steady  state  responses  of  large  electrical  networks.  Of  particular 
interest  were  the  delay  times  for  circuit  activation,  power  consumption 
during  the  trajectory,  and  detonator  firing  at  the  termination  of  the 
trajectory. 

Fuze  circuit  operation  was  investigated  by  considering  projectile  firings 
from  selected  altitudes  and  effective  muzzle  velocities.  The  thermoelectric 
open-circuit  voltage  (as  a function  of  elapsed  time  from  firing)  for  these 
different  conditions  was  used,  and  a total  of  eight  situations  were  analyzed. 


* Developed  at  Kirtland  Air  Force  Base. 
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The  revised  circuit  shown  in  Figure  4 was  modeled  for  SCEPTRE  analysis.  An 
arbitrary  number  of  300  thermocouples  was  selected  for  the  thermoelectric 
voltage  source.  The  calculated  thermoelectric  voltage  as  a function  of  time 
after  firing  used  in  these  analyses  is  given  in  Figure  17.  The  calctu. 
lations  were  terminated  at  2.5  seconds  after  firing,  except  for  the  low  altitude 
case.  The  firing  angle  was  maintained  constant  at  -20  degrees  from  the  hori- 
zontal axis  for  all  cases  considered. 

Table  IV  illustrates  the  relationships  among  the  firing  parameters  and  the 
operation  of  the  proximity  fuze  circuit.  The  time  delay  column  indicates 
when  the  circuit  becomes  active  after  firing.  The  duration  column  indicates 
the  elapsed  time  after  firing  that  the  circuit  is  actively  sensing  a target. 

In  other  words,  the  proximity  fuze  is  no  longer  active  if  the  thermoelectric 
voltage  has  dropped  below  the  threshold  voltage  of  approximately  9 volts 
necessary  to  operate  the  RF  circuitry.  However,  the  detonator  circuitry 
remains  armed  after  the  RF  sensor  ceases  to  function  due  to  the  stored  energy 
in  the  2pP  capacitor.  The  detonator  is  then  activated  upon  impact  through 
the  closing  of  the  impact  switch,  SI. 

From  Table  IV  it  can  be  seen  that  the  revised  proximity  fuze  circuit  is 
"compatible"  with  the  thermoelectric  voltage  source.  From  the  use  of 
220  couples,  the  length  of  time  the  proximity  fuze  circuit  remains  active 
is  more  than  2.5  seconds  at  muzzle  velocities  of  3500  ft/sec  or  more.  This 
duration  corresponds  to  an  active  range  for  the  proximity  fuze  of  more  than 
6000  feet. 

To  extend  the  time  of  operation  (hence,  range)  for  the  proximity  fuze  cir- 
cuitry, the  following  steps  may  be  taken: 

• Increase  projectile's  initial  velocity 

• Increase  number  of  thermocouples 

• Utilize  energy  storage 

• Revise  the  circuit  design. 

All  or  any  of  these  factors  can  be  used  to  increase  the  active  life  of  the 
circuit.  From  Table  IV  we  can  see  that  in  most  cases  the  delay  time  (time 
until  activation)  is  satisfactory  for  operational  conditions.  However, 
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the  duration  of  an  active  proximity  circuit  can  be  improved  by  increasing 
the  effective  muzzle  velocity.  For  example,  comparing  the  thren  initial 
velocity  conditions  at  a 3000-foot  altitude  indicates  that  an  initial  velocity 
of  3700  ft/sec  or  greater  provides  a longer  duration,  hence,  range. 

Alternatively,  the  number  of  thermocouples  could  be  increased  to  366  (maximum 
for  space  available),  which  would  provide  additional  voltage.  The  voltage 
would  increase  by  a factor  of  1.66. 

Energy  storage  was  also  considered  in  extending  the  active  life  of  the 
circuit.  Initial  calculations  indicate  a capacitance  on  the  order  of  1000  uF 
would  be  needed  to  supply  the  necessary  current  for  1 second.  However,  the 
space  limitations  associated  with  the  present  circuit  design  prevent  the 
use  of  chip  capacitors  on  the  surface  of  the  existing  circuit  boards.  As 
an  alternative,  a separate  substrate  containing  a multiplic  ity  of  chip 
capacitors  or  a discrete  element  capacitor  might  be  added  to  the  projectile 
tip. 

Also,  the  RF  circuit  may  be  further  modified.  The  circuit  was  not  optimized 
since  the  objective  was  to  determine  if  mating  a proximity  fuze  circuit  with 
a thermoelectric  power  supply  was  feasible.  With  the  construction  of  a 
breadboard  circuit  and  further  computer  circuit  analyses,  the  proximity  fuze 
circuit  may  be  optimized  to  function  for  the  different  firing  condi- 
tions using  less  than  the  present  input  power  at  40  milliwatts. 

8.  ADVANTAGES  AND  LIMITATIONS  OF  ENERGY 
STORAGE  IN  RF  PROXIMITY  FUZES 

The  previous  discussion  has  indicated  that  one  method  for  extending  the  period 
during  which  the  RF  fuze  is  "active"  is  to  use  energy  storage  By  using 
storage,  excess  power  produced  by  the  thermoelectric  converter  during  the 
early  phase  of  the  trajectory  would  be  available  to  supply  the  needed  thresh- 
hold  voltage  and  current  during  the  later  phase  of  the  tra jectory . . (For 
example,  see  Figure  18,3000-ft  altitude,  4000-ft/sec  muzzle  velocity.) 
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Based  on  the  findings  of  this  study,  the  following  are  requisite  conditions 
for  the  effective  use  of  energy  storage  with  20-ms  projectiles: 

* Aircraft  altitude  (during  firing)  should  be  below 
25,000  ft. 

a Effective  muzzle  velocity  should  be  at  least 
3500  ft/sec. 

Several  approaches  to  energy  storage  can  be  considered — thermal  and  electrical 
energy  storage.  Thermal  storage  is  effected  through  the  use  of  thicker  hot 
straps  (see  Figure  15)  or  surrounding  the  thermoelectric  converter  with  an 
outer  layer  of  metal  or  ceramic.  For  example,  a 0.010  in.  thick  sleeve  of 
metal  (copper  or  aluminum)  enclosing  the  disc  modules  (see  Figure  7)  would 
affect  the  output  power  versus  time  as  shown  in  Figure  19. 

A second  possible  approach  to  energy  storage  involves  the  storage  of  elec- 
trical energy.  As  indicated  in  Subsection  7,  a 1000  pF  capacitor  would  be 
required  to  extend  the  "active  life"  of  the  RF  fuze  by  1.0  second,  or  about 
500  pF  of  capacitance  for  0.5  second  increase  in  active  life.  "Chip" 
capacitors  which  can  be  placed  in  planar  substrates  (consistent  with 
the  present  circuit  configuration)  are  available  in  capacitances  up  to  only 
2 pF.  Hence,  250  to  500  of  these  0.225  in.  x 0.210  in.  chips  would  be 
required  for  the  present  projectile.  Alternatively,  a discrete  tubular 

O 

capacitor  component  may  be  used  instead  but  will  occupy  a volume  of  0.5  in. 
to  1.0  in.^  for  the  500  to  1000  pF  capacitors,  respectively.  For  the 
case  of  the  20-mm  projectile,  this  volume  requirement  limits  the  use  of 
electrical  storage  unless  (1)  the  projectile  tip  can  be  enlarged  or  (2)  the 
central  core  of  the  projectile  is  bored  out  to  accommodate  a tubular  capacitor. 
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• Altitude  = 3000  ft 


FIGURE  19.  TEMPERATURE  DIFFERENCE  ACROSS  TOE  THERMOELECTRIC 
CONVERTER  FOR  SELECTED  HOT  STRAP  THICKNESSES 


SECTION  VI 


STRUCTURAL  ANALYSIS  OF  DISC  THERMOELECTRIC  CONVERTOR 

A preliminary  structural  analysis  of  the  projectile  tip  as  modified  for  a thermo- 
couple generator  is  discussed  next.  Since  the  detailed  design  of  the  projectile 
will  not  be  initiated  until  Phase  II,  these  efforts  have  considered  only  the 
overall  structural  integrity  of  the  thermopile  with  an  attempt  to  identify  possible 
problem  areas. 

1 . ANALYTICAL  APPROACH 

The  overall  dimensions  and  configuration  of  the  tip  are  shown  in  Figure  7.  Load- 
ings to  the  tip  are  a spin  ratio  of  2000  rev/sec  and  a peak  acceleration  level 
of  120,000  g's.  In  addition,  aerodynamic  heating  gives  a (222C)  temperature  rise 
at  the  surface  of  the  Vespel  polyimide  discs  which  make  up  the  thermocouple 
generator.  It  was  assumed  that  the  impulse  accelerating  the  projectile  is  such 
that  dynamic  stress  wave  propagation  effects  are  negligible.  The  120,000  g 
acceleration  can  then  be  considered  a quasistatic  inertia  loading. 

It  is  understood  that  the  projectile  prior  to  the  present  modification  has  been 
demonstrated  to  withstand  the  imposed  loadings.  Thus  it  is  necessary  only  to 
consider  the  region  in  the  area  of  the  circumferential  groove  in  which  the 
thermoelectric  convertor  is  positioned.  The  following  properties  were  used  for 
the  polyimide  substrate  in  the  structural  analysis: 

Elastic  Modules  * 4.8  x 10^  lb^in.2, 

Density  = 0.046  lb/in. 3, 

Strength  ■*  10,000  psi, 

Thermal  Expansion  Coefficient  ■ 2.85  x 10"  5 in./in./F.  s K 

It  was  recognized  that  the  polyimide  material  has  an  elastic  modulus  about  50 
times  less  than  that  of  the  adjacent  steel  components.  Thus  it  was  reasonable  to 
assume  that  loads  from  the  projectile  tip  will  be  transferred  entirely  through  the 
steel  shell,  and  not  through  the  plastic  discs.  Also  it  was  assumed  that  the  steel 
shell  can  be  treated  as  rigid  relative  to  the  polyimide. 
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2. 


ROTATIONAL  STRESSES 


Although  the  discs  will  probably  be  bonded  or  mechanically  fixed  into  the  groove 
of  the  steel  body,  in  the  stress  analysis  no  credit  was  taken  for  support  of  the 
disc  by  the  body.  Hius  rotational  stresses  were  calculated  by  treating  the  discs 
as  freely  spinning  rings.  The  maximum  hoop  stress  in  such  a ring  is  given  by 

(Ts)max  = (~£-)ow  (b  ), 

where 

b = outside  radius  = 0.3  inch 
a = inside  radius  = 0.2  inch 

3 

c = density  = 0.046  lb/in. 
w = speed  of  rotation  = 12,566  rad/sec 
v = Poissan's  ratio  = 0.3 

Evaluating  the  above  equation  for  these  parameters  gives  the  tensile  hoop  stress 
as 

(Js)max  = 15  28  lb/in. 

2 

which  is  much  less  than  the  10,000  lb/in.  tensile  strength  of  the  polyimide  material. 
The  strain  associated  with  this  stress  is 

E = J/E  = (1528/4.8  x lO"*)  x 100  percent 
= 0.32  percent. 

The  bismuth  telluride  layer  would  be  subjected  to  essentially  this  level  of  strain. 
Iliis  is,  however,  an  upper  limit  on  strain  level,  since  no  support  of  t lie  discs 
by  the  steel  shell  has  been  considered.  In  fact,  the  imposed  strain  should  be 
much  less  than  this  value. 

3.  ACCELERATION  STRESSES 

The  stresses  due  to  the  120,000  g acceleration  were  estimated  by  considering  the 
weight  of  the  0.3  inch  high  stack  of  discs,  and  distributing  this  force  over  t he 
annular  cross  sectional  area  of  these  discs.  The  weight  was  determined  as  0.00216 
lb.  At  120,000  g 1 8 this  becomes  260  lb.,  and  when  distributed  over  the  base  of 
the  discs  gives  a stress  of  1655  psi.  This  is  much  less  than  the  disc  strength 
of  about  10,000  psi. 
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4. 


THERMAL  STRESSES 


A upper  bound  on  the  thermal  stress  in  the  polyimide  discs  can  be  obtained  from 
the  equation 

J - E xoCT  ■=  (4.8  x 105)  (2.85  x 10'5)  (400) 

* 5472  psi , 

which  is  less  than  the  reported  strength  of  the  material.  Perhaps  of  greater 
concern  is  possible  cracking  of  the  bismuth  telluride  layer  due  to  differential 
thermal  expansion  relative  to  the  polyimide  substrate.  The  substrate  expansion 
coefficient  of  2.85  x 10  ^ in./in./F  is  relatively  high  with  respect  to  bismuth 
telluride,  which  is  on  the  order  of  1.0  x 10"^  in./in./F.  Unless  the  expansion 

of  the  polyimide  discs  is  restrained  by  the  steel  shell,  the  bismuth  telluride 
layer  will  be  placed  in  tensile  stress. 

Unfortunately,  efforts  to  obtain  mechanical  properties  of  bismuth  telluride  were 
unsuccessful.  Lacking  such  data,  one  would  be  unable  to  determine  whether 
bismuth  telluride  can  withstand  a 400  F (222C)  temperature  cycle  when  subject  to 
strains  induced  by  a polyimide  substrate. 

5.  COMMENTS  ON  STRUCTURAL  INTEGRITY 

Details  of  the  design  and  assembly  of  the  thermoelectric  generator  device  have 
not  yet  been  formulated.  It  is  recommended  that  a more  nearly  complete  stress 
study  be  made  of  the  final  configuration  (to  be  developed  in  Phase  II  of  this  study). 
The  preliminarv  structural  analysis  indicates  that  the  polyimide  substrate  will 
no'  be  overstressed  for  the  anticipated  service. 
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SECTION  VII 


DYNAMIC  STABILITY  ANALYSIS 


The  stability  of  a spin-stabilized  projectile  is  a function  of  its  spin 
angular  momentum,  the  aerodynamic  forces  and  moments,  and  the  mass  properties. 
The  spin  angular  momentum  is  determined  by  the  spin  rate  and  the  axial  moment 
of  inertia.  The  aerodynamic  forces  and  moments  are  explicitly  related  to  the 
exterior  shape  of  the  projectile,  Mach  number,  dynamic  pressure,  Reyno.' ds 
number,  and  location  of  the  center  of  gravity.  Thus  it  is  seen  that  the  mass 
properties  influence  both  the  angular  momentum  and  the  aerodynamics.  In  this 
brief  analysis  the  effect  of  mass  property  variations  on  the  stability  of  a 
20mm  proximity  fuze  projectile  are  to  be  considered.  Unfortunately,  these 
variations  are  related  in  a complex  manner  to  the  other  factors  influencing 
stability  mentioned  above. 


In  part 'cular.  this  analysis  is  to  consider  the  effects  on  stability  of  four  mass 
properties,  viz.,  axial  moment  of  inertia  (1^),  transverse  moment  of  inertia 
(I),  mass  (m) , and  distance  from  the  nose  to  the  center  of  gravity  (X  ). 

y ^8 

There  are  two  ways  of  viewing  the  stability  and  the  influence  of  mass  property 
changes.  In  the  first,  the  effect  on  the  damping  coefficients  Ap  and  \n  (see 
page  35  of  Reference  2)  is  considered.  Stability  is  attained  when  both 

coefficients  are  negative.  For  small  mass  property  increments  it  is  sufficient 

to  investigate  the  coefficients  with  the  smallest  magnitude.  For  the 

cases  to  be  considered  below  only  the  effect  on  >.p  is  studied.  The  other 

stability  criterion  utilizes  a plot  of  the  generalized  stability  factor  (s  ) 

d 

and  the  gyroscopic  stability  factor  (s  ).  The  plot  is  usually  given  as  s 

8 8 

and  8^  for  the  abscissa,  and  ordinate,  respectively.  By  definition  the  pro- 
jectile is  stable  if  the  points  on  the  plot  are  contained  within  an  envelope 
defined  by  the  equation 

1 


-1 


8 


sd^2_Sd^' 
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These  two  means  for  viewing  stability  are  equivalent.  The  advantage  of  the 
first  is  that  it  gives  quantitative  information  on  the  distance  required  for 
a disturbance  to  dampen  to  some  fraction  of  its  initial  amplitude.  The  second 
stability  criterion  can  give  insight  into  whether  the  projectile  is  near  the 
boundary  of  the  stability  regime. 

Instead  of  determining  the  stability  for  one  or  two  sets  of  mass  property 

parameters,  a more  general  approach  will  be  taken.  vaiue  0f  \ % 8 an(j 

P 8 

s^  will  be  given  at  sea  level  conditions  for  two  Mach  numbers  for  the  original 
proximity  fuze  missile  designed  by  Motorola.  In  addition  the  partial  derivatives 

of  these  three  quantities  will  be  shown  with  respect  to  1^,  1^,  M and  X /L 
(where  L is  body  length).  This  will  allow  one  to  obtain  the  effect  on  stability 
of  an  arbitrary  but  small  variation  in  the  mass  properties. 

The  reason  for  considering  only  sea  level  conditions  is  that  stability  problems 

will  arise  there  first.  Since  dynamic  pressure  will  decrease  with  altitude 

(if  speed  is  constant),  the  projectile's  stability  will  increase.  This  is 

because  the  sources  of  instability  are  certain  ones  of  the  aerodynamic  moments  (with 

static  pitching  moments  being  the  major  one),  and  these  moments  are  proportional  to  the 

dynamic  pressure. 

Table  V contains  a list  of  the  partial  derivatives  to  be  used  in  determining 
X p at  Mach  numbers  1 and  4.  A new  value  of  Xp  is  determined  by  the  formula 

4 

Xp  " kp„omln.l  * 2 Af J 

where  ^nominal  = va^ue  damping  coefficient  for  Motorola  RF 

fuze  design  at  appropriate  Mach  number,  per 
caliber  (see  Table  5) 

I 8 lug- ft2 

Iy  SlUg-ft2 

m slugs 

X /L  -nondimens ional 
eg 

^new  "^nominal. 
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TABLE  V. 

PARTIAL  DERIVATIVES  FOR  DETERMINING  Xp  (SEA  LEVEL) 


Mach 

No. 

1 

4 


Nominal  Values  iXp  , rt4  fcXp  , „4  &Xp  , „4  fcXp  ,^4 
of  Xp  jjf  * 10  Sf  X 10  V X 10  $0?  7L)  X 10 

-0.757  x 10“  0.135  x 106  0.00432  x 106  -113.6  -0.385 

-0.471  x 10"4  -0.795  x 106  0.0134  x 106  251  -5.6 


TABLE  VI. 

PARTIAL  DERIVATIVES  FOR  DETERMINING  S "l  (SEA  LEVEL) 


Mach 

No. 


1 

4 


Nominal 

Value 

0.0792 

0.635 


)S 

T? 


-1 


%■ 


fcs  -1 

n&r 


-0.0543  x 10" 
-0.221  x 106 


0.00698  x 10 
0.0121  x 106 


0 

0 


^S  1 
.ft 
i(X  /L) 
£ft 

-0.3438 

-1.19 


TABLE  VII. 

PARTIAL  DERIVATIVES  FOR  DETERMINING  S,  (SEA  LEVEL) 

a 


Mach  Nominal  ^d  ^Sd  ^Sd  ^Sd 

No.  Value  Si”  ST”  TnT  $(X  71) 

* 1 ££ 

1 0.702  -0.0261  X 106  0.00902  x 106  0 -1.896 

4 0.634  -0.104  x 106  0.0070  x 106  -61.1  -1.45 


TABLE  VIII. 

EFFECT  OF  MASS  PROPERTY  CHANGES  ON  THE  STABILITY  PARAMETERS 


Mach  No. 

Original  Values 

x s -i 

S, 

X 

New  Values 

s -1 

S . 

P g 

d 

p 

K 

d 

1 

_4 

-0.757  x 10  0.0792 

mm 

-0.757  x 

10"4  0.0802 

0.715 

4 

-0.471  x 10-4  0.635 

u 

-0,466  x 

10-4  0.603 

0.629 
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The  nominal  values  and  partial  derivatives  required  to  calculate  new  values  of 

Sg-1  and  S . are  found  in  Tables  VI  and  VII.  The  formula  given  above  is  used  with 

either  S or  S.  in  place  of  X . 

g d P 

Figure  21  shows  a plot  of  versus  and  the  points  for  the  Motorola  RF 

fuze  design  at  Mach  numbers  1 and  4. 


For  the  new  projectile  with  the  thermoelectric  powered  fuze,  the  following  mass 
property  increments  were  calculated 

AT  = 0.08x10"”^  slug-Ft^ 

AI  = 0.51x10“^  slug-Ft^ 

= 0.129x10”^  slug 

AX  /L  = -0.00532 
cs 

where  a positive  increment  denotes  that  the  new  projectile's  value  of 

that  parameter  is  higher  than  the  one  with  the  battery.  The  old  and  new 

values  of  X_,  S “*■,  and  S.  at  Mach  numbers  1 and  4 are  given  in  Table  VIII.  It 
■ 8 ^ 

is  seen  that  the  effect  of  the  design  modifications  is  practically  negligible. 
In  fact  no  attempt  was  made  to  plot  the  new  values  in  Figure  20  since  the  points 
would  be  almost  on  top  of  the  original  ones. 
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SECTION  VIII 


COST  EFFECTIVENESS  OF  THERMOELECTRIC  CONVERTER 

The  purpose  01  this  task  was  to  explore  alternative  thermopile  manufacturing 
approaches.  In  view  of  the  size  of  the  thermoelectric  elements  required  for 
this  application,  thin  film  fabrication  techniques  will  be  required.  Several 
techniques  for  vapor  deposition  of  bismuth  telluride  including  flash  vapor 
deposition  and  DC  sputtering,  were  considered.  The  DC  sputtering  technique  was 
selected  for  use  in  preliminary  studies  based  on  the  higher  degree  of  control 
afforded  and  the  potential  for  mass  production  manufacturing  techniques. 

The  principal  disadvantage  of  sputtering  is  that  it  requires  6 to  12  hours 
to  deposit  0. 002-in. -thick  thermoelectric  elements.  However,  commercial 
sputtering  equipment  would  permit  depositing  bismuth  telluride  on  up  to 
10,000  disc  modules  simultaneously. 

The  steps  involved  in  the  evaluation  of  the  cost  effectiveness  of  thin  film 
manufacturing  methods  included  investigating  alternative 

• Substrate  (disc)  materials 

• "Masking"  techniques 

• Hot-strap  and  cold-strap  deposition  methods 

• Sputtering  parameters. 

Vespel  was  identified  as  a suitable  substrate  material  for  use  in  the  fabrica- 
tion of  disc  modules  during  Phase  II  of  this  study.  Vespel  was  selected  because 
of  (1)  its  ability  to  wi^^tand  exposures  t-o  temperatures  of  up  to  500  F,(260C), 
(2)  its  availability  ir  *■  .e  required  thickness  (0.035  in.),  and  (3)  its  low 
thermal  conductivity. 

Studies  also  revealed  that  the  dimensions  of  the  conductive  straps  (which 
connect  the  p-  and  n-type  thermoelectric  elements  and  wrap  around  the  edge 
of  the  substrate)  are  large  enough  that  it  may  be  possible  to  deposit  thj 
straps  by  the  standard  silk-screen-and-fuze  process.  Otherwise,  thermal 
evaporation  of  conductive  metal  through  a mask  may  have  to  be  employed  to 
deposit  the  straps. 

The  investigation  of  both  sputtering  parameters  and  the  effectiveness  of 
alternative  masking  techniques  was  undertaken  and  involved  preliminary 
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sputtering  trials  First,  bismuth  telluride  bar  stock  (supplied  by  Melcor) 
was  diced  into  1 in.  x 2 in.  x 0.062  in.  sections  in  order  to  make  up  a 4 in. 
diameter  "target"  for  use  in  BCL's  sputtering  equipment.  This  bismuth  telluride 
target  provided  for  the  deposition  of  0.001  in.  to  0.002  in.  thick  layers  of 
bismuth  telluride  on  selected  substrates.  These  preliminary  deposition  studies 
provided  an  important  input  to  the  study  of  cost  effective  thermopile  manu- 
facturing approaches  as  well  as  direction  for  the  prototype  converter 
development  scheduled  for  Phase  II  of  this  study. 

Next,  Bi^Te^  (p-type)  was  sputtered  onto  glass  slides  by  using  stainless  steel 
masks.  Two  identical  masks  were  made  with  the  EDM  (Electrical  Discharge 
Machine)  method.  The  overall  size  of  the  masks  was  1 x 2 x 0.008  in.  The 
slots  cut  in  the  mask  was  0.100  in.  long  x 0.025  in.  wide.  Three  of  these 
slots  had  a space  of  d.065  in.  between  them  and  the  fourth  had  a space  of 
0.020  in.  Another  form  of  mask  which  consisted  of  metal  strips  0.25  in.  wide  x 
1.0  in.  long  and  0.003  in.  thick  was  used.  These  were  clamped  over 
glass  slides  with  a space  between  the  strips  of  0.125  in. 

The  parameters  for  sputtering  the  B^Te^  were  as  follows: 

Cathode-Substrate  distance  — 2.5  in. 

Cathode  material  — 1 in.  x 1 in.  Bi^Te^  (P-tyPe) 
positioned  on  top  of  a 4-in. -diameter  iron  plate 
Base  pressure  — 5 x 10  ^ torr 
Air  pressure  — 40  p, 

Foreline  pressure  — 200  p, 

Cathode  voltage—  2.5  Kv 
Cathode  current  — 40  ma 
Power  — 100  watts 

2 

Cathode  power  density  — 8 watts/in. 

Deposition  rate  — 280  X/rain 
Deposition  time  — 15  hr. 

Measurements  of  the  coating  thickness  were  made  with  a Talysurf.  As  was  expected, 
the  coating  through  the  mesk  was  rounded  near  the  edges. 
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Sputtered  bismuth  telluride  elements  of  dimensions  and  separation  distance 
representative  of  the  actual  disc  converter  are  illustrated  in  Figure  21. 

As  indicated  above,  a stainless  steel  "physical"  mask  was  used  in  the 
sputtering  of  these  thin  film  thermoelements.  The  significance  of  these 
results  is  that  less  costly  physical  masking  techniques  can  be  used  in  the 
thin-film  deposition  of  bismuth  telluride. 

Electrical  property  measurements  were  made  in  order  to  characterize  the 
sputtered  thin  film  thermoelements.  These  measurements  were  made  on  0.125 
in.  wide  x 0.80  in.  long  films  which  ranged  in  thickness  from  0.0005  in. 
to  0.001  in.  The  results  of  preliminary  measurements  indicate  that  the 
Seebeck  coefficient  of  sputtered  films  is  within  15  percent  of  that  of  the  bismuth 
telluride  source  material.  The  electrical  resistivity  measurements,  however, 
indicate  that  the  films  have  a factor  of  two  or  more  higher  elec- 
trical resistivity  than  the  bulk  material.  Microscopic  inspection  of  the  thin 
films  indicated  that  most  of  the  samples  exhibited  microcracks.  This  micro- 
cracking likely  contributed  to  the  high  electrical  resistivity.  Analysis 
of  the  sputtering  system  revealed  that  the  glass  slides  probably  heated  up 
to  200  C (380  F)  or  more  during  the  sputtering  process.  This  high  tempera- 
ture combined  with  the  different  thermal  expansion  rates  of  the  bismuth 
telluride  and  the  glass  slide  likely  led  to  excessive  stresses,  and,  hence, 
cracking  in  the  thermoelectric  material.  The  only  elements  that  were  not 
cracked  were  the  smaller  0.25  in.  wide  x 0.100  in.  long  elements.  However, 
these  smaller  elements  were  deposited  in  a second  sputtering  series  which 
did  net  provide  representative  p-type  bismuth  telluride  films.  This  was 
due  to  the  fact  that  the  "target"  was  fabricated  by  using  bismuth  telluride 
wafer  segments  which  were  not  in  good  thermal  contact  with  the  target  base 
and,  hence,  the  water  cooled  heat  sink.  By  overheating,  the  more  volatile 
species  of  bismuth  telluride  evaporated  during  the  earlier  sputtering  trial-- 
thereby  disproportioning  the  bismuth  telluride  alloy.  This  problem  can  be 
readily  overcome  by  fabricating  a monolithic  bismuth  telluride  target  and 
achieving  good  therms'  contact  with  the  water  cooled  heat  sink. 
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13X 

FIGURE  21.  PHOTOGRAPH  OF  SPUTTERED  BISMUTH  TELLURIDE  FILMS  (O.lOO  IN 
LONG  X 0.025  IN  WIDE  X 0.0005  IN  THICK) 


A preliminary  estimate  has  been  made  of  the  costs  of  fabricating  the  selected 
disc  thermoelectric  modules  Involving  bismuth  tellurlde.  These  very  preliminary 
estimates  are  sunmarlzed  In  Table  9 and  Indicate  that  the  total  thermoelectric 
convertor  may  be  fabricated  for  less  than  $1.80  per  convertor  (containing  7 
modules).  These  estimates  have  been  made  conservative  by  assuming  a twc  year 
payback  on  capital  equipment,  200  working  days  per  year,  two  runs  per  day  per 
set  of  equipment  where  , where  a run  equals  10,000  disc  convertors.  Further  refinement 
of  the  cost  analysis  is  not  warranted  due  to  the  uncertainty  in  ultimate 
deposition  processes  and  substrate  material.  However,  it  does  indicate  an 
approximate  figure  for  the  convertor  concept  in  its  present  stage  of  development. 
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TABLE  IX. 

APPROXIMATE  COST  BREAKDOWN 
FOR  DISC  THERMOELECTRIC  CONVERTORS 


Materials 

Labor 

Unit  Disc 

Item 

Cost , $ 

Cost, $ 

Cost,$ 

Bismuth  Telluride 

200/lb. 

— 

0.0266 

Vespel  Substrate 

10/in.3 

— 

0.0600 

Sputtering  Equipment* 

— 

— 

0.0375 

Electrode  Deposition 
Equipment 

— 

— 

0.012 

Supporting  Labor** 

— 

$15/hr. 

0.108 

Building  Rental, 
Other  Costsx 

Total 

Cost  Per 

Convertor 

$ 1.79 

* Assume  $3000,000  for  equipment  plus  automatic  mask/disc 
handling  equipment 

+ Assume  $100,000  for  thick  film  equipment. 

**  Three  workers  per  shift,  three  shifts, 
x Based  on  $100,000  per  annum  costs. 
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SECTION  IX 


SAFETY  ANALYSIS 

The  aerodynamical ly  heated  thermoelectric  convertor  concept  developed  in  this 

study  has  also  been  analyzed  from  the  standpoint  of  safety  during  ground 

handling,  storage,  transportation,  and  firing.  One  existing  feature  of  the 

projectile  is  the  "spin-back"  safety  measure,  which  prevents  arming  the 

projectile  until  the  rotational  effects  of  firing  "arm"  the  firing  pin.  In 

addition  to  this  feature,  the  aerodynamically  heated  thermoelectric  power 

supply  offers  an  additional  measure  of  safety  to  its  users.  This  second  safety 

feature  derives  from  the  fact  that  the  RF  proximity  fuze  is  not  activated  until 

a threshold  voltage  (and  current)  is  developed  within  thermoelectric  convertor. 

In  order  to  develop  the  necessary  voltage/current  output  from  the  convertor,  a 

2 2 

sufficiently  large  heat  flux  (>10.0  Btu/sec.  ft  or  13.0  watts/cm  ) must  be 
applied  and  sustained  for  at  least  50  to  100  msec. 

In  the  event  the  thermoelectric-powered  RF  proximity  fuzes  were  exposed  to  a 

fire  (e.g.,  often  simulated  by  a 1450  F (790  C)  radiant  source),  the  heat 

2 

flux  would  reach  only  about  6 watts/cm  . This  heat  flux  is  only  about  one 
half  the  heat  flux  required  to  develop  the  threshold  voltage  necessary  to 
activate  the  RF  proximity  fuze.  Hence,  in  the  event  of  a fire,  the  projectile 
tip  will  heat  up  without  developing  a significant  temperature  gradient  within 
the  thermoelectric  converter.  At  some  point,  a fire  will  heat  the  RF  proximity 
fuze  components  to  temperatures  beyond  the  operational  limits  of  the  electronic 
components.  For  example,  the  Vespel  substrates  will  thermally  decompose  above 
500  F (260  C) . The  electronic  components  in  the  RF  circuit  will  fail  rapidly 
as  the  temperatures  approach  these  temperatures. 

Another  measure  of  safety  provided  by  the  thermoelectric-powered  fuze  is  the 
safe  arming  distance.  This  distance  is  dictated  by  the  thermal  inertia  of  the 
thermoelectric  convertor  in  the  presence  of  the  imposed  aerodynamic  heating 
rate.  Based  on  analyses  discussed  in  Part  7 of  Section  VI  (See  Table  4),  the 
safe  arming  interval  is  at  least  50  msec  or  about  150  feet. 
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The  incidence  of  "hung  rounde"  will  not  activate  the  thermoelectric-powered 
RF  proximity  fuze  unless  the  temperature  of  the  projectile  environment  reaches 
or  exceeds  1800  F (990  C) . Heating  or  cooling  of  the  projectile  prior  to 
firing  will  not  significantly  affect  the  time  history  of  the  temperature 
difference  developed  across  the  convertor  via  aerodynamic  heating.  However, 
the  ultimate  hot  junction  of  the  convertor  will  depend  on  the  initial 
projectile  temperature.  It  should  be  noted  that  the  output  voltage  developed 
has  a first  order  dependence  on  the  temperature  difference  across  the 
convertor.  Therefore,  since  the  safe  arming  distance  depends  on  the  time 
history  of  the  convertor  output  voltage,  it  follows  that  the  safe  arming 
distance  will  be  relatively  insensitive  to  the  initial  temperature  of  the 
projectile. 
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SECTION  X 


CONCLUSIONS  AND  RECOMMENDATIONS 


During  Che  seven  month  period  of  this  Phase  I study,  a number  of  analyses  have 

been  performed  aimed  at  assessing  the  feasibility  of  utilizing  aerodynamical ly 

(2) 

heated  thermoelectric  convertors  to  power  previously  developed'  ’ RF  proximity 
fuzes.  The  collective  results  of  this  study  indicate  that  such  a thermoelectric 
power  supply  is  feasible  and  can  be  made  compatible  with  the  existing  RF 
proximity  fuze — in  terms  of  both  the  fuze  voltage/current  requirements  and 
the  safe  arming  distance  desired. 

The  principal  findings  from  this  study  include: 

• 20  mm  RF  proximity  fuze  may  be  powered  with  as  little  as 
40  mw  at  an  input  voltage  of  7.5  volts. 

• Bismuth  telluride  is  the  most  attractive  thermoelectric 
material  for  the  number  of  couples  required  and  the 
threshold  aerodynamic  heat  flux  required. 

• Disc  thermoelectric  module  concept  may  offer  a cost 
effective  means  for  distributing  200  to  400  couples  over 
the  available  surface  area  on  the  20-mm  projectile. 

• Aircraft  should  have  an  airspeed  of  at  least  80  mph  and 
not  more  than  760  mph  in  order  to  provide  2.5  seconds  of 
active  power  to  the  RF  fuze  circuit  while  flying  at 
3000  feet  and  firing  at  -20  degrees. 

• At  significantly  higher  altitudes  (>25,000  feet),  the 
above  airspeed  limits  will  be  somewhat  lower  but 
trajectory  llfetimer  of  up  to  3.5  seconds  may  be  possible — 
without  any  energy  storage. 

• At  increasingly  lower  altitudes  and  airspeeds  above  80  mph, 
there  is  an  increasing  incentive  to  use  energy  storage  to 
provide  excess  energy  produced  during  early  phase  of  trajectory 
for  later  phase  of  trajectory. 
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• Electrical  (capacitor)  storage  Is  one  means  for  extending  the 
active  life  of  the  RF  proximity  fuze.  However,  physical  size 
limitations  within  20-mm  projectile  tip  provide  increases  in 
trajectory  of  only  0.5  seconds. 

• Preliminary  structural  analysis  indicates  that  the  polyimide 
substrate  will  not  be  overstressed  for  the  anticipated 
service.  However,  the  effects  of  differential  thermal 
expansion  between  the  Vespel  and  the  bismuth  telluride 
thermoelectric  elements  pose  some  uncertainties  at  this  time. 

• Dynamic  stability  analyses  of  the  modified  projectile  tip 
indicate  that  the  changes  in  center  of  mass,  moment  of  inertia, 
and  mass  of  projectile  tip  are  sufficiently  small  to  have  a 
negligible  effect  on  the  dynamic  stability  of  the  overall 
projectile. 

• Preliminary  safety  analysis  indicates  that  the  thermoelectric 
convertor  requires  a sufficiently  high  heat  flux  (to  produce 
threshold  power  for  fuze)  to  minimize  the  chance  that 
accidental  fires  or  "hung"  rounds  will  activate  fuze. 

• Preliminary  thin  film  deposition  studies  indicate  that 
bismuth  telluride  can  be  sputtered  in  the  required 
dimensions  by  using  relatively  low  cost  physical  masking 
techniques.  Electrical  property  measurements  on  these 
preliminary  trial  depositions  revealed  the  Seebeck  coefficient 
to  be  within  15  to  20  percent  of  that  of  the  bulk  material. 

The  electrical  resistivity,  however,  was  a factor  of  two  or 
more  larger  than  that  of  the  bulk  material.  The  high 
resistivity  of  these  samples  has  been  traced,  in  part,  to 
microcracks  in  the  thin  film  specimens — possibly  due  to 
differential  thermal  expansion  of  bismuth  telluride  with  respect 
to  the  glass  substrate. 

• 30  mm  RF  proximity  fuze  can  be  powered  by  aeroheated  thermo- 
electric convertors  for  trajectory  durations  of  at  least 
3.5  seconds  owing  to  its  larger  projectile  mass. 
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Based  on  the  above  findings,  it  is  concluded  that  aerodynamically  heated 
thermoelectric  convertors  provide  a feasible  source  of  power  for  RF  proximity 
fuzes  in  20  mm  antj  30  mm  projectiles.  Certain  technical  uncertainties  exist, 
however,  at  this  phase  of  the  concept  development: 

• Which  deposition  method  is  best  suited  to  this  task  upon 
considering  the  need  for  high  volume/low  cost  production 
and  high  figure  of-merit  material? 

• Can  bismuth  telluride  in  0.1  in.  lengths  withstand  the 
differential  thermal  expansion  effects  associated  with 
the  use  of  Vespel  substrates? 

• Can  acceptably  low  resistivity  material  be  deposited 
without  the  need  for  annealing  (sintering)  following 
deposition? 

• What  effects  will  long  term  storage  of  20  mm  projectiles 
have  on  the  thermoelectric  convertors?  Will  moisture 
and  heat  and  cold  have  any  effects? 

On  the  basis  of  the  findings  of  this  Phase  I study,  it  is  recommended  that  a 
Phase  II  study  be  initiated  emphasizing 

• Refinement  of  thin  film  deposition  of  bismuth  telluride 
on  polyimide  substrates  including  enhanced  cooling  of 
target  and  substrate  during  sputtering* 

• Characterization  of  deposited  bismuth  telluride  thermo- 
elements, before  and  after  thermal  cycles  simulating 

2 

several  seconds  at  a heat  flux  of  13  to  15  watts/cm  . 

• Optimization  of  deposition  parameters,  substrate 
preparation,  and  electrode  configuration. 

• Evaluation  of  prototype  projectile  tips  incorporating  disc 
thermoelectric  modules  including  output  voltage,  current, 
and  hot  and  cold  Junction  temperatures  as  functions  of 
time  during  1.5  to  2.5  second  simulated  trajectories. 
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APPENDIX  A 

DESCRIPTION  OF  THE  COMPUTER  PROGRAM  TRUMP 
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APPENDIX 


DESCRIPTION  OF  THE  COMPUTER  PROGRAM  TRUMP 

TRUMP  solves  a general  nonlinear  parabolic  partial  differential 
equation  describing  flow  in  various  kinds  of  potential  fields,  such  as  fields 
of  temperature,  pressure,  and  electricity  and  magnetism;  simultaneously,  it 
will  solve  two  additional  equations  representing,  in  thermal  problems,  heat 
production  by  decomposition  of  two  reactants  having  rate  constants  with  a 
general  Arrhenius  temperature  dependence.  Steady  state  and  transient  flows 
in  one,  two,  or  three  dimensions  are  considered  in  geometrical  configurations 
having  simple  or  complex  shapes  and  structures.  Problem  parameters  may 
vary  with  spatial  positions,  time,  or  primary  dependent  variables — temperature, 
pressure,  or  field  strength.  Initial  conditions  may  vary  with  spatial 
position,  and  among  the  criteria  that  may  be  specified  for  ending  a problem 
are  upper  and  lower  limits  on  the  size  of  the  primary  dependent  variable, 
upper  limits  on  the  problem  time  or  on  the  number  of  time  steps  or  on  the 
computer  time,  and  attainment  of  steady  state.  Solutions  may  be  obtained 
by  use  of  explicit  or  implicit  difference  equations,  or  by  an  optimized 
combination  of  them. 
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